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ABSTfU\CT 


llio  influence  of  combined  states  o(  stress  on  the  shear  response 
.ilonp,  material  piincipal  directions  in  off-axis,  unidirect  iiinal  ,  composite 
cotipons  is  examined  for  systems  wliose  matrix  material  obeys  the  von  Mises 
yield  condition.  Such  analysis  is  motivated  by  trends  in  experimental 
data  observed  for  at  least  two  composite  systems  that  indicate  deviation 
trom  pure  shear  behavior  alonp.  material  directions  for  various  off-axis 
con f  i>;u rat  ions . 

T\\o  yield  function  for  plane  stress  of  a  transversely  isotropic 
compiisite  lamina  consist  in>;  of  stiff,  linearly  elastic  fibers  and  a 
von  Mises  rruitrix  material  is  formulated  in  terms  of  Hill’s  elastic  stress 
concent  rat iiMi  factors  and  a  single  plastic  constraint  parameter.  The  above 
are  sul>sequent ly  evaluated  on  the  basis  of  observed  average  lamina  and 
constituent  resptnise  for  the  Avco  5505  horon-epoxy  system.  It  is  shown 
that  inclusion  of  residual  stresses  in  Lite  yield  function  together  with  the 
incorporat ion  of  Pubey  and  Hill  let’s  concept  of  generalized  yield  stress 
lor  anisotn'pic  media  in  the  constitutive  equation  correctly  predicts  the 
trends  observed  in  exjter  iments .  Tlie  incorporat  ion  of  the  strong  axial 
stress  interaction  necessary  to  predict  the  correct  trends  in  the  shear 
respi'use  is  directly  traced  to  the  high  residual  axial  stresses  in  the 
matrix  induced  during  fabrication  of  the  composite. 
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(50/.  fiber  volume  fraction)  under  combined  loads. 
Normal: shear  load  ratio  is  8:3,  ref.  3 


Experimental  sliear  stress-strain  response  along 
materialprincipal  directions  for  various  off-axis 
tension  coupons.  AVCO  5505  boron/epoxv  system  ... 


Predicted  shear  stress-strain  response  along  material 
principal  directions  for  various  off-axis  tension 
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Predicted  shear  stress-strain  response  along  rotated  axes 
for  an  equivalent  von  Mises  isotropic  coupon  in  tension. 
Based  on  pure  5505  boron/epoxy  system  . . 


Yield  shear  stress  of  AVCO  boron/epoxy  system  along  the 
material  principal  directions  of  an  off-axis  tension 
coupon  as  a  function  of  the  fiber  orientation  ... 


redicted  shear  stress-strain  response  along  material 
principal  directions  for  various  off-axis  tension  coupons 
based  on  Dubey  and  Hlllier*s  generalized  yield  stress 
concept  and  micromechanical  considerations  including 
residual  stresses.  AVCO  5505  boron/epoxy  system  . 
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INTRODUCTION 


Th6  use  of  the  off~axis  tensile  specimen  for  determination  and/or 
verification  of  elastic  properties  and  fracture  strength  of  unidirectional 
c» mposltc  materials  is  well  documented^  eq.  Sih  [Ij*  As  most  composite 
materials  exhibit  varying  degrees  of  nonlinear  response  prior  to  failure 
in  an  off-axis  tension  test,  this  specimen  would  also  appear  well  suited 
for  verification  of  the  various  assumptions  in  a  number  of  nonlinear 
theories  recently  proposed  in  the  literature,  as  well  as  some  of  their 
consequences.  An  indirect  outgrowth  of  certain  assumptions  employed  in 
one  such  theory,  for  instance,  is  a  recently  proposed  method  for  determi¬ 
nation  of  nonlinear  lamina  shear  stress-strain  response  on  the  basis  of 
the  10®  off-axis  tension  test  12]. 

The  nonlinear  response  of  unidirectional  composites  is  a  complex 
phenomenon  which  can  be  caused  by  a  number  of  diverse  mechanisms  such  as 
nonlinear  constituent  behavior,  damage  accumulation  through  fiber  or 
matrix  cracking  or  debonding,  fiber  rotation  or  any  combination  of  the 
above  to  mention  the  better  known  causes.  The  extension  of  results  of 
the  uniaxial  test  along  principal  material  directions  to  constitutive 
models  for  multiaxial  loading  situations  is  not  always  straightforward 
due  to  the  possibility  of  various  stress  interactions  taking  place  at  the 
microlevel.  Foy  [3]  for  example  demonstrated  the  Influence  of  the  trans¬ 
verse  normal  stress  component  on  the  shear  strain  response  (and  vice 
versa)  along  the  material  axep  for  composites  with  stiff,  linearly  elastic 
fibers  and  ductile  matrices  using  a  finite  element,  micro-mechanical  anal¬ 
ysis.  He  showed  that  for  combined  normal  shear  loading  in  a  fixed  ratio 
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yielding  initiated  aooner  and  the  extent  of  inelastic  behavior  was  tiore 
pronounced  than  for  either  normal  or  shear  stress  acting  alone  (Fig.  1). 

A  number  of  macromechanical  constitutive  models  have  been  proposed  in 
the  literature  [4-8]  to  describe  the  nonlinear  response  of  lamina  with 
stiff,  elastic  fibers  and  significantly  more  compliant  and  often  ductile 
matrices.  The  early  attempts  [4,5]  disregarded  the  possibility  of  stress 
interaction  in  the  nonlinear  range  during  combined  loadings,  while  in  the 
later  formulations  this  effect  was  taken  into  account,  often  in  an 
heuristic  manner.  One  such  model  due  to  Hashin  et  al  [6]  is  reminiscent 
of  the  total  deformation  theory  of  plasticity  in  that  the  stress  interac¬ 
tion  is  accounted  for  through  a  loading  function,  based  upon  an  Invariant 
formulation,  which  is  used  to  determine  the  different  functional  forms  of 
the  various  strain  components.  The  loading  function  is  limited  to  quad¬ 
ratic  terms  in  transverse  normal  and  shear  stresses  :hus  eliminating  the 
possibility  of  the  axial  normal  stress  influence  on  the  transverse  and 
shear  strain  behavior. 

Data  generated  by  Cole  and  Pipes  [9]  are  employed  in  this  paper  to 
re-examine  and  study  the  extent  of  stress  interaction  in  off-axis  unidi¬ 
rectional  boron-epoxy  coupons  with  the  aid  of  Kashin's  model  and  a  micro¬ 
mechanical  approach  utilizing  a  yield  function  concept  based  on  Hill’s 
method  [10]  for  determining  the  elastic  stress  concentration  factors. 

The  results  of  the  micromechanical  analysis  are  subsequently  employed  in 
Dubey  and  Hillier’s  generalized  yield  stress  formulation  [11]  to  study 
the  nonlinear  shear  stress-strain  response  along  the  principal  material 
directions  in  off-axis  tension  coupons.  Several  fiber  orientation  are 
considered  including  the  10°  coupon  which  has  been  proposed  for  the 
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STRESS-STRAIN  CURVES  FOR  UNiniRECTIONAI.  BORON/EPOXY  (50Z  FIBER 
VOLUME  FRACTION)  UNDER  COMBINED  LOADS.  NORMAL: SHEAR  LOAD  RATIO 
IS  8:3,  REF.  3. 
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der.ermlnation  of  pure  sh^ar  response. 

Th.  lanu»n«  of  .cros.  Intor.ctlon  on  tbo  -boor  veapono.  1»  tbo 

p„0  abear  ra.ponao  obfafoad  fao™  aa.dfaaaf fonal  aoba,.  ^aso  aap.rf.aaf.f 
abov  tbaf  tb.  aafabf  of  In.U.tfc  babav.of  la  .of  a  -.ofo.la 
,„„cflon  of  flbaa  ofia.faao..  Tba  l.alaatlc  aba.f  .trains  (tor  a  spacl- 
fled  abaar  afraa.)  daaraaaa  1.  «S.lr.da  aa  tba  flb.r  a.sla  la  l.craaaad 
fro.  15-  to  30-.  but  aa  tba  flbar  a.gla  la  l.ar.aaad  furtbar  to  45  tba 
inalaatla  .train.  Inaraaaa.  (Sl.Uar  babavlor  naa  baan  obaarvad  b,  tba 
antbora  for  tba  pr.pblta-polpl.lda  off-anla  conponaO  ^a  of  tb.  .oa.a 
tbla  papar  1.  to  provlda  an  axpl.n.tlon  for  aucb  babavlor. 
in  tba  first  part  of  tba  papar  tba  ..dal  of  Haabln  at  .1  la  brlafly 
ontllnad.  It  la  tban  naad  to  pradlct  tba  abaar  atraaa-atraln  r.apona, 

.loop  .atarlal  araa  nalng  nnldlractl.nal  data  t.kan  fro.  Raf.  191.  *« 

a„plrlaa.  ob.arvatlon  la  ..da  tb.t  .»dlflcatlon  of  H.abln’a  .od.l  by 
incorporation  of  an  axial  nor.al  atr.as  Intaractlon  tar.  significantly 
Lprovaa  corralatlon  nltb  tba  axparl.ant.l  data.  ...ad  on  tba  abov. 
raanlta.  „lcro.acb.nlcal  fo.-.l.tlon  of  tba  yield  function  for  tba 
cpoalta  in  tba  p.-sanca  of  stiff  fiber.  1.  e.Ployad  along  vltb  bubaV. 
generalized  yield  stress  concept  to  examine  the  effect  of  t.ie 
component  on  yielding  in  the  matrix  phase  and  its  consequences  on  the 
„onlinear  shear  response  along  material  directions.  It  is  shown  that 
incorporation  of  tba  prasance  of  rasldnal  curing  stresses  In  t 
leads  to  tb.  correct  prediction  of  trend.  In  tbe  exparLantal  .bear  atre.a- 

Strain  data  of  Cole  and  Pipes. 


STRESS  INTERACTION  MODEL 


Fo  lowltiB  t'l;issical  concepts  of  deformation  theory  of  plasticity 
Hashln  ec  al  [6J  assume  that  the  inelastic  part  of  strain  can  be  expressed 
in  the  following  manner, 

summation)  (1) 

where  L  is  some  general  quadratic  loading  function  of  the  stresses. 

Since  composites  which  employ  stiff  fibers  such  as  boron  or  graphite 
behave  linearly  (or  nearly  so)  under  loads  applied  in  the  fiber  direction, 
the  plastic  response  is  assumed  to  be  ] ^mited  to  transverse  and  shear 
strains.  Furthermore,  the  effect  of  the  stress  in  the  fiber  direction 
on  the  remaining  strain  components  is  neglected  and  thus  Eq.  (1)  is 
expressible  as  follows: 


12  2G^2  ^12  °21  ^®12^ 

P  ^22  2  2 

22  “  E^  ^22  ®22  ^°12^ 


(2) 


where  the  interaction  term  Oj^2®22  excluded  on  the  basis  of 

material  symmetry  arguments  and  only  strains  relevant  to  laminate  analysis 
retained.  The  functional  form  of  ^22^^^  determined  from 

uniaxial  L*  nsion  and  torsion  tests  along  material  directions.  It  is  noted 
that  no  explicit  mention  is  made  of  the  yield  function  and  its  relationship 
to  the  loading  function  in  this  formulation  despite  apparent  separation  of 
the  total  strain  into  elastic  and  inelastic  or  plastic  portions.  Combining 
plastic  and  elastic  portions  of  strain,  the  relevant  total  strain  components 


become : 
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where  the  constants  M,  N,  a 


and  T  are  obtained  from  the  Ramberg-Osgood 


approximation  of  uniaxial  tests 


THEORETICAL-- EXPERIMENTAL-CORRELATION 


Data  presented  by  Cole  and  Pipes  on  the  Avco  5505  boron-epoxy  system 
have  been  employed  to  determine  the  Ramberg-Osgood  parameters  in  Hashin’s 
model.  These  are  given  in. the  appendix.  Predicted  shear  response  in  the 
material  coordinate  system  for  the  various  off-axis  fiber  orientations  is 
plotted  in  Fig.  3.  While  Hashin’s  model  predicts  monotonically  decreasing 
response  with  increasing  off-axis  angle,  the  experimental  data  points  to  a 
reversal  around  30®  off-axis  angle  (Fig.  2> . 

The  concept  of  stress  interaction  in  an  elastoplastic  material  can  be 
illustrated  by  considering  the  shear  stress-strain  response  along  axes 
Inclined  at  some  angle  with  respect  to  the  loading  direction  of  an  isotropic 
tensile  specimen.  The  hypothetical  isotropic  material  is  taken  to  have  the 
same  shear  properties  as  the  considered  boron-epoxy  system  for  illustration 
purposes.  Using  the  isotropic  theory  (in  this  case  the  flow  and  defor¬ 
mation  approaches  are  coincident)  the  shear  response  along  the  inclined 
axes  is  given  by: 

N-1 

^  ^12  ^  12  f  2  1  .  2.  ^  2.,,  2  ,12, N-1 

h2  2^*  2077  f  3  +  3  ® 

12  J  2  y 

with  yield  governed  by: 

°12  ^  f  T  cotan^e)]  =  S^,  O<0<9O*  (5) 

/3 

(or  alternatively,  “2  ^  sin20)  where  S  is  the  yield  strecs  in 

in  pure  shear.  Eq.  (4),  plotted  in  Fig.  4,  predicts  monotonically 
increasing  shear  response  along  rotated  axes  with  increasing  angle  between 
0®  and  45®  with  reversal  occuring  at  45®.  This  follows  from  the  symmetry 
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ENGR.  SHEAR  STRAIN 


PREDICTED  SHEAR  STRESS-STRAIN  RESPONSE  ALONG  MATERIAL 
PRINCIPAL  DIRECTIONS  FOR  VARIOUS  OFF-AXIS  TENSION 
COUPONS.  HASHIN  ET  AL  MODEL.  AVCO  5505  BORON-EPOXY 
SYS! 


SHEAR  STRESS  (KSI) 
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FIG.  4.  PREDICTED  SHEAR  STRESS-STRAIN  RESPONSE  ALONG 
ROTATED  AXES  FOR  AN  EQUIVALENT  VON  MISES 
ISOTROPIC  COUPON  IN  TENSION.  BASED  ON  PURE 
5505  BORON-EPOXY  SYSTEM. 


II 


of  the  yield  function  (Eq.  5)  about  the  45*  angle.  For  comparison  purposes 
Hashin's  theory  (Fig.  3)  yields  the  following  relation  for  the  shear  res¬ 
ponse  along  material  axes  Inclined  at  the  same  angle  as  the  reference  axes 
in  the  hypothetical  Isotropic  specimen 


’12 


'12 


12  2G 


12 


2G 


12 


*  9 

(1  +  (-^)  tan  el 


N-1 

2 


y 


(6) 


The  absence  of  the  resolved  shear  stress  reversal  in  the  above  model  c.-i  be 
directly  traced  to  the  absence  of  the  axial  stress 


Interaction  term. 


MICROMECHANICAL  CONSIDERATIONS 


Assuming  that  the  matrix  can  be  modelled  as  a  von  Mises  solid  that 
controls  the  nonlinear  response  of  the  composite,  the  yield  function  of 
the  lamina  can  be  determined  from  the  knowledge  of  the  stress  concentration 
factors  at  yield.  We  take  the  view  that  these  stress  concentration  factors 
are  related  to  the  average  elastic  concentration  factors  through  a  single 
plastic  constraint  parameter.  This  appears  to  be  partially  substantiated 
by  various  finite-element,  micromechanical  studies,  c.f.  Dvorak  and  Rao 
[12],  which  have  revealed  that  although  yielding  usually  start 
flber/™.trlx  lntbrf»c..  tl»  .tr.s.  v.btor  M.rly  rbdl.1  under 

proportional  loading.  The  yield  zone  on  the  other  hand  spreads  rapidly 
throughout  the  entire  .atrlx  »lth  Increasing  external  tractions. 

n.e  elastic  stress  concentration  factors  can  be  determined  fro.  the 
knonledge  of  the  average  la.lna  as  veil  as  constituent  tesponse  foUovlng 
Hill's  outline  110).  The  relationship  betu.en  average  natrlx  and  external 
stresses  at  yield.  Including  residual  stresses,  1.  then  assumed  to  be 
expressibl6  as  follows: 


^llm 


^22in 


^12m 


=  K 


^11 

^21 

0 


B 


B 


12 

‘22 

0 


0 

0 


B 


‘66  — * 


(7) 


and  the  elastic 


.  iTd.  ...  m'lfrix  r  residual  and  no  subscript  external 

where  subscript  m  refeis  to  matrix,  r 

^  '  o, 

stresses.  Also,  “  °33mr  ®22mr  “  ®33mr 

stress  concentration  factors  . are  given  in  the  appendix. 

Here,  is  the  previously  mentioned  plastic  constraint  factor  to  be 
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determined  from  experimental  data,  and  the  single  scale  factor  A  is  a 
consequence  of  transverse  isotropy  of  the  composite. 

Employing  the  above  formulation,  von  Mises  yield  function  for  the 
transversely  isotropic  lamina  in  plane  stress  becomes: 


3  3  ^®12“®12®22''‘®22^”22‘^p 

+  3  (1  -  A^Oiin,r^2®ir®21^°ll’^p  +  3  “  A^‘^llrar^^®12“®22^°22’^p 


(8) 


where  K  is  the  yield  shear  stress  of  the  matrix  material. 

The  importance  of  the  residual  stresses  cannot  be  overemphasized. 
Employing  a  single  inclusion  composite  cylinder  model  and  temperature 
dependent  material  properties  for  this  particular  system,  the  average  cur¬ 
ing  stresses  have  been  estimated  to  be: 


Although  time  dependent  response  has  not  been  considered  due  to  insuf¬ 
ficient  data,  the  reduction  in  residual  stresses  due  to  viscoelastic 
response  will  be  partially  off-set  by  the  use  of  a  single  inclusion  com¬ 
posite  cylinder  model  which  generally  underestimates  induced  stresses. 

The  above  yield  function  (Eq,  8)  is  plotted  in  Fig.  5  in  terms  of  the 
resolved  shear  stress  along  the  material  axes  and  the  off-axis  angle. 


NORMALIZED  SHEAR  YIELD  STRESS 
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Also  plotted  are  the  yield  functions  of  an  equivalent  isotropic  material 
and  a  transversely  isotropic  lamina  with  no  residual  stresses*  It  is 
clearly  seen  that  the  reversal  of  the  shear  stress-strain  response  evident 
in  Coleys  and  Pipes*  data  can- be  explained  on  the  basis  of  initial  yielding 
in  the  matrix.  To  do  this  we  formulate  a  generalized  yield  stress  function 
according  to  Dubey*s  and  Hillier*s  suggestion  as  follows: 

and  assume  that  the  plastic  shear  strain  is  a  function  of  o.  The  coeffi¬ 
cients  C..,  D...  j  are  determined  from  the  yield  function  given  by  Eq.  (8) 
and  the  pure  shear  response*  These  are  given  in  the  appendix.  The  result¬ 
ing  shear  stress-strain  curves  along  material  axes  for  the  various  off-axis 
angles  are  presented  in  Fig*  6  indicating  the  extent  of  stress  interaction. 
It  is  seen  that  the  use  of  the  10®  off-axis  specimen  for  determination  of 
pure  shear  response  can  lead  to  noticeable  deviations  for  a  material  system 
that  can  be  modelled  by  the  assumed  set  of  governing  equations* 


SHEAR  STRESS 


ENGR.  SHEAR  STRAIN 

Fir,  6  PREDICTED  SHEAR  STRESS-STRAIN  RESPONSE  ALONG  MATERIAL  PRINCIPAL 
DIRECTiSs  FOR  VARIOUS  OFF-AXIS  TENSION  COUPONS  BASED  ON  DUBEY 
AND  HILLIER’S  GENERALIZED  YIELD  STRESS  CONCEPT  AND  MICRO¬ 
MECHANICAL  CONSIDERATIONS  INCLUDING  RESIDUAL  STRESSES.  AVCO 
5505  BORON-EPOXY  SYSTEM. 


CONCLUSIONS 


The  deviations  from  pure  sliear  response  along  material  axes  in  an  off- 

axes  tension  test  for  the  Avco  5505  boron-epoxy  system  have  been  explained 

to  a  large  extent  on  the  basis  of  micromechanical  formulation  of  a  yield 

surface  and  a  generalized  yield  stress  concept.  The  model  is  applicable 

to  systems  with  stiff,  linearly  clastic  fibers  and  ductile  matrices  whose 

response  can  be  approximated  by  that  of  a  von  Mises  solid.  The  trends 

evident  in  the  data  of  Cole  and  Pipes  for  the  shear  response  along  the 

mater-^al  principal  directions  are  traced  to  the  strong  axial  stress  (o..) 

^  ii 

interaction  which  in  turn  follows  from  relatively  significant  residual 
stresses.  The  model  predicts  that  in  the  absence  of  these  stresses  the 
Influence  of  axial  stress  on  nonlinear  shear  stress-strain  response  for 
low  o£f-axis  fiber  orientations  is  only  significant  in  the  0®-5®  off-axis 
range.  When  the  residual  stresses  are  included,  however,  this  influence 
is  noticeable  up  to  approximately  15®.  This  is  interesting  in  view  of  the 
recently  proposed  test  method  for  determination  of  pure  nonlinear  shear 
stress-strain  response  on  basis  of  the  10®  off-axis  tension  test. 
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APPENDIX 


i)  Ramberg-Osgood  parameters  for  Che  Avco  5503  Boron-Epoxy  system 


M 

Oy(psi) 

N 

Ty(p8i) 

3.A34 

26,155.2 

4.430 

9,593.9 

ii)  Elastic  stress  concentration  factors 
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where 


G  *  0.19  X  10  psi 
m  ^ 

E  =  0.49  X  10^  psi 
m  ^ 

V  =  0.31 
m 

=  0.50 


*  24.0  X  10  psi 
=  58.0  X  10^  psi 
■=  0.20 


^12  *  ^  IG  psi 

=  30.1  X  10^  psi 


"22 


=  2.87  X  lo'*  psi 


Vj2  =  0.225 

iii)  Determination  of  generalized  yield  stress.  From  experiment 
.P 
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66  ^12  contracted  notation.  Thus 
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From  longitudinal  tension  and  compression  we  obtain:  ^  In 

contracted  notation,  where 


?  - 


‘■w 

(x«p 


the  magnitudes  of  yield  stress  in  tension  and  compression 
respectively.  Similarly  in  transverse  tension  and  compression 


C2  =  n 


where 


<Y«^) 


L 


The  parameters  S,  and  n  have  been  evaluated  from  the  yield 

function  given  by  Eq.  (8)  and  the  experimental  data.  Tlie  plastic 

constraint  factor  K  has  been  determined  on  the  basis  of  yielding  (i.e. 
P 

deviation  from  linearity)  in  transverse  tension. 


